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n^ . Abstract. We calculate nuclear suppression J^^iyi of heavy quarks produced from the 

initial fusion of partons in nucleus-nucleus collisions at RHIC and LHC energies. We 
take the shadowing as well as the energy loss suffered by them while passing through 
Quark Gluon Plasma into account. We obtain results for charm and bottom quarks at 

O ' several rapidities using different mechanisms for energy loss, to see if we can distinguish 

2 ' between them. 

> 

00 ' 1- Introduction 

o' 

The heavy ion colhsion experiments at RHIC and LHC are designed with a hope to 
ly^ . explore the existence of a new form of matter known as Quark Gluon Plasma (QGP) and 

^ . to explore its properties. The estimation for the energy density [H |2l [3l H] attained in 

these collisions using the Bjorken formula [5] is well beyond the energy densities where 
QGP is expected to be formed. The temperatures reached at RHIC, as revealed from 
K^ , several studies (see e.g., O [TJ |8] for a compilation) are also much larger than the values 

provided by Lattice QCD calculations for the critical temperature for a transition to 
QGP [9j. Strong confirmation of the formation of the QGP is given by observation of 
a large elliptic flow [TO], jet-quenching |Tl], and the recombination of partons as the 
mechanism of production of hadrons at intermediate transverse momenta [12]. Still 
higher temperatures are likely to be reached at LHC. 

Heavy quarks are produced from the initial fusion of gluons {gg — )■ QQ) or light 
quarks (gg — )■ QQ). This pair would be produced at r ^ 1/2Mq ^ 0.1 fm/c. Their 
large mass ensures that their production can be treated using pQCD and that it is 
nearly negligible at later times. These will traverse the QGP, colliding with quarks and 
gluons and radiating gluons before appearing as charm or bottom mesons or baryons. 
Thus the final spectra for these hadrons would contain information about the energy 
loss suffered by the heavy quarks. The unique importance of heavy quarks as probes of 
QGP lies in their large mass. This leads to a considerably reduced production of heavy 
quarks in comparison to light quarks and gluons which are produced copiously. The 
strong interaction during the collision conserves flavour. Even at the LHC, the reverse 
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process of annihilation of heavy quarks {QQ — )■ gg, etc.) can be safely ignored. Thus 
heavy quarks and in turn charm and bottom hadrons will stand out in the back-ground 
of a multitude of light hadrons, and one can in principle track them. 

Along with other reasons, an interest in the study of energy loss of heavy 
quarks was triggered by the large back-ground that correlated charm or bottom decay 
provides [131 IH] to the thermal dileptons which have been considered a signature of the 
formation of QGP for a long time [IH [161 Il3 [HI US] ■ It was pointed out by Shuryak [20] , 
Lin et al. [211 ES], Kampfer et al. [23], and Mustafa et al. [2l| 125] that the correlated 
charm and bottom decay could be suppressed if the energy loss suffered by heavy quarks 
before they form D or B mesons was accounted for. Since then several attempts have 
been made to estimate the energy loss of heavy quarks as they proceed through QGP. 
Of course the possibility to identify the vertex of D or B meson decay will further enrich 
this study. 

These results have been put to a rigorous test by the measurement of single electrons 
from heavy ion collisions at RHIC which show a clear evidence for the loss of energy 
by the heavy quarks [26]. Their possible thermalization is also indicated by the elliptic 
flow that they show [27] . 

As indicated earlier, the temperature likely to be reached at LHC in collision of 
heavy nuclei could be even larger and thus this energy loss will play a more significant 
role. The opening of a much wider window in rapidity at LHC is also likely to provide 
widely differing media at different rapidities through which the heavy quarks would 
propagate. 

Thus a valuable test of various theories for energy loss suffered by heavy quarks 
can be performed by studying it at RHIC and LHC and at different rapidities. 

We study these effects in terms of nuclear modification factor Raa for heavy quarks. 
In these initial studies we calculate the average energy loss suffered by them as they 
pass through the QGP and compare the resulting pt distribution with the same for 
proton-proton (pp) collisions to get Raa- Since the mass of charm or bottom quarks is 
quite large, the pt distribution of these quarks will closely reflect the pt distribution of 
D or B mesons. 

We employ a local fluid approximation [13l [HI [28] in order to picture the medium 
at larger rapidities. We shall come back to this later. 

The paper is organized as follows. As we need to compare the spectra of the heavy 
quarks from relativistic heavy ion collisions with those for pp collisions, as a first step 
we study the heavy quark production in LO pQCD and compare our results with a 
NLO pQCD calculation. We find that single quark distribution calculated using LO 
pQCD supplemented with a K-factor adequately reproduces the NLO results as well as 
the available experimental data. Next we estimate the average energy loss suffered by 
heavy quarks of a given energy using various mechanisms discussed in the literature. 

Finally we perform a Monte Carlo calculation to obtain the average change in the 
transverse momentum spectra of heavy quarks for nucleus-nucleus collisions and get 
Raa as a function of px for different rapidities. We add that this work is not intended 
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as a complete review and the readers may see Ref . [2^ [20] , for other treatments. 

More detailed calculations where the consequences of energy loss of heavy quarks on 
the correlated charm or bottom decay and modification of the back-to-back correlation 
of heavy quarks are discussed, will be published shortly. 

2. Heavy quark production in pp collisions 

At lowest order in pQCD, heavy quarks in pp collisions are produced by fusion 
of gluons {gg — )■ QQ) or light quarks (gg — )■ QQ) [31]. The so-called flavour excitation 
process {qQ -^ qQ and gQ — )■ gQ) is now known to be suppressed when the NLO 
processes are taken into account [32], [33l [3l] . In addition, Brodsky et al. [35l [36] have 
shown that the total contribution of intrinsic charm in the midrapidity region is small 
even though most of the heavy quarks are produced in this region. 

The cross-section for the production of heavy quarks from pp collisions at lowest 
order is given by 



2xiX2PT^[fP {xi,Q'^)ff\x2,Q'^)crij{s,i,u) + 
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ff\x,,Q^)fP{x2,Q^)a,,{sXu)]^^ + ^^3)■ (1) 



In the above equation, i and j are the interacting partons, // and /j are the partonic 
structure functions and Xi and X2 are the fractional momenta of the interacting hadrons 
carried by the partons i and j. The relation between pt and fractional momentum xi 
or X2 through their respective rapidities can be written as 

xi = !^ (e^i + e^^), X2='^ {e-y^ + e"^^), (2) 

yS ys 



where m^ is the transverse mass, J M"^ + p"^, of the produced heavy quark. The 
function a = da/dt, the short range subprocess for the heavy quark production is 
defined as: 

|A^P for the heavy quark production processes gg^QQ and qq^QQ are expressed 
through the mass of the heavy quark and Mandelstam variables s, t, and u as 
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The running coupling constant Us at lowest order is 

127r 



a,9 



(5) 



(6) 
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where A^/ = 3 is the number of active flavours and A = Aqcd- We use the factorization 
and renormalization scales as Q'^ = m'^. We refer the readers to Vogt et al. [38j for 
results on variations of these scales. We also carry out the calculation of differential 
cross section for heavy quarks in pp collision at NLO in pQCD using the treatment 
developed by Mangano, Nason, and Ridolfi (MNR-NLO) [39j. All the calculations are 
carried out by neglecting the intrinsic transverse momentum of the partons. 

The effect of nuclear shadowing in high energy nucleus-nucleus collisions is well 
known [lOl SH |l2l HS]. With the increase of the mass number of the nucleus 
and increasing contribution of terms having small x, the effect becomes more 
pronounced. We introduce the shadowing effect in our calculations by using EKS 98 
parameterization [Hj for nucleon structure functions. We take CTEQ4M [15] structure 
function set for nucleons. 

We shall see that x dependence of the shadowing function introduces interesting 
structures in the nuclear modification factor as a function of pt, y, and the incident 
energy, because of the large mass of the quarks. 

In Fig. [1] we compare our results for heavy quark px distribution obtained using 
lowest order pQCD for pp collision with the results from NLO-MNR calculation 
at midrapidity for charm and bottom quarks at RHIC and LHC energies. These 
comparisons suggest a K factor of ~ 1.5 - 3 for our lowest order calculations for agreement 
with NLO results. 

In Fig. [2] we compare the bottom quark production cross-section obtained using 
lowest order pQCD in pp collision at 630 and 1800 GeV energies with UAl, CDF and 
D0 data [l6l HTj. We find a good description of these data using lowest order pQCD 
with a K-factor. 

We calculate the total cross-section for charm quark production at lowest order for 
the process pp— )-cc as a function of y/s considering the charm quark mass as 1.2 GeV 
and 1.6 GeV. We have also included results for M^ (3 GeV) = 0.986 GeV [IH], suggested 
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Figure 1. [Left panel] Comparison of our lowest order pQCD results with the NLO- 
MNR calculation for charm quark [Mc = 1.6 GeV) at midrapidity. [Right panel] Same 
for bottom quark [Mb — 4.7 GeV) at midrapidity. 
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Figure 2. Comparison of our result for a (pT > p™*") for production of bottom 
quarks with experimental data. 



recently. Though the last value is not at the pole Q = Mc, it may serve as the lower limit 
to the mass of the charm quark. In the left panel of Fig. [3] we compare these results and 
the results from NLO-MNR calculation with the experimental data points ^7\ HH |50] . 
Here also our lowest order pQCD calculation show a good agreement with experimental 
data points for Mc = 1.2 GeV. In the right panel of Fig. Owe compare our results with 
the NLO-MNR calculations up to ^5 = 15000 GeV. 

We also calculate the total cross-section for bottom quark production for pp— > bb 
as a function of y/s considering the bottom quark mass as 4.7 GeV and 4.163 GeV |18]. 
In Fig. Hlwe compare these results and the results obtained from NLO-MNR calculation 
with the experimental data points [51]. This comparison is quite impressive as at 
Mf, = 4.7 GeV our lowest order result accurately reproduce the result of NLO-MNR 
calculation with K = 2.5. 
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Figure 3. Total cross-section for pp— >-cc compared witli experimental data at varying 
Vs. 



Thus we see that the pt distribution and production cross-section for the heavy 
quarks calculated in lowest order pQCD and supplemented with a K-factor ^ 2.5 
reproduces the results at NLO for pp collisions. In view of this we feel that these 
distribution would be adequate for calculating Raa- 
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Figure 4. Total cross-section for pp— > bb compared witti experimental data varying 



3. The Initial conditions and the Evolution of the Plasma 



The heavy quarks produced at the initial stage pass through the QGP, where they loose 
energy by colliding with quarks and gluons and also by radiating gluons. The energy loss 
will depend upon the path-length of the heavy quarks in the plasma, the temperature 
evolution of the plasma, and the energy and mass of the heavy quarks. 
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In order to proceed we make several simplifying assumptions. It is expected that 
the heavy quarks will loose most of their energy when the temperature is still large, 
i.e. during the earliest times after the formation of QGP. During these early times, we 
can neglect the transverse expansion of the plasma. We assume a Gaussian rapidity 
density distribution for the particles produced and further assume that the initial 
rapidity distribution of the quarks and gluons follows this distribution closely. This 
would correspond to an isentropic expansion of the plasma at all rapidities, and should 
be sufficient for our initial studies. 

Thus we assume the rapidity distribution of the density of gluons as [131 El]- 

We take (^) ?^900 and a = 3 for Au+Au collisions at RHIC |52] and ^3300 and 

^ dy 'q ! i 

a = 4 for Pb+Pb colhsions at LHC [53]. 

The Bjorken cooling is then assumed to work locally at different rapidities, and 
we consider the passage of a heavy quark having rapidity y in a fluid having an 
identical fluid rapidity. This approximation, which corresponds to assuming a boost- 
invariant expansion along with a local fluid approximation, has been used earlier 
in literature [131 El I2H]- A more complete study would use a (3 + 1) dimensional 
hydrodynamics [SH |55] , which we plan to use in future publications. 




Figure 5. The distance, L, covered by a heavy quark while passing through the QGP. 
For central collisions the results for (L) will not depend on $. 



We consider a heavy quark produced in a central collision, at the point (r, $), and 
moving at an angle with respect to f in the transverse plane. In general the distance 
covered by the heavy quark before it exists the QGP, will vary from to 2R, where R 
is the radius of the colliding nuclei. The distance covered by the heavy quark in the 
plasma, L, is given by [56j: 



R2 — r^ sin^ — r cos 0. (8) 
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We can estimate the average distance travelled by the heavy quarks in the plasma as: 



{L) 



/rdr fL(0, r)TA^(r,6 = 0)d, 


/rdr fT^^(r,6 = O)d0 





(9) 



In the above the nuclear overlap function T^aI^, b = 0) provides the probability of 
production of heavy quarks in hard binary collisions. We find that (L) is 5.78 fm for 
Au+Au collisions at RHIC and 6.14 fm for Pb+Pb collisions at LHC, and is about 20% 
smaller than the radii of the colliding nuclei, as the appearance of the nuclear overlap 
function gives a larger weight to the points having smaller r. 

As the heavy quarks loose most of their energy in interaction with gluons, it is 
enough to consider only the distribution of gluons. Their density at the time r can be 
written as [57] : 



1 dNn 
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The corresponding temperature [57], assuming a chemically equilibrated plasma is 
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The rapidity dependence of the temperature of the plasma at a typical r = (L)/2 at 
RHIC and LHC is given in Fig. |6l 
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Figure 6. Variation of temperature of QGP with rapidity y at a typical r. 



Assuming that the QGP is formed at ro = 0.2 fm/c, we estimate the Tq at y = 
for RHIC as 377 MeV and LHC as 555 MeV. More detailed studies do suggest a larger 
formation time of ~ 0.5 fm/c, which will correspondingly reduce the initial temperatures. 
This will not affect our results as we approximate the expanding and cooling plasma. 
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with one at a temperature determined at r = Lcfj/2, which is much larger (see the 
following discussion). Further assuming, Bjorken's cooling law, T^r = constant, this 
provides that the plasma would cool down to the transition or critical temperature 
Tc~ 160 MeV by rc~ 2.6 fm/c at RHIC and 8 fm/c at LHC One can easily calculate 
the corresponding values at larger rapidities. 

Considering the velocity of the quark as vt=Pt/^tT'T, it would take a time tl = 
(L)/vt to cross the plasma. If Tc > t^ the heavy quark would be inside QGP during 
the entire period, tq to r^. However, if Tc < tl, only while covering the distance 
vt X Tc, would the heavy quark be in the QGP phase. We further approximate the 
expanding and cooling plasma with one at a temperature of T at r = (L)eflf/2, where 
(L)cff = min [(L), vt x r^]. This procedure has been used frequently |l57j . 

4. Mechanisms for Energy Loss 

Next we discuss the energy loss mechanisms that we have included. As mentioned 
above, repeatedly, the heavy quarks loose energy both by collisions as well as radiation 
of gluons. A number of formalisms have been proposed for the collisional as well as the 
radiative energy loss of heavy quarks in the literature. We shall consider the following 
treatments for the collisional energy loss. 

Bjorken [58] has considered the collisional energy loss of light quarks as analogous 
to the energy loss of a charged particle passing through a medium and losing energy 
by ionizing the medium. His expression for massless quarks was adapted by Braaten 
and Thoma to the case of heavy quarks [see Eq. lA.l] . We shall continue to label this 
mechanism as Bjorken for clarity. Braaten and Thoma (BT) [591 125] also modified the 
expression for the energy loss suffered by muons while traversing QED plasma, to obtain 
the collisional energy loss of a heavy quark as it passes through the QGP [see Eqs. IA.2I 
and IA.3l in the Appendix A]. These results are valid for collisions where the momentum 
transfer q << E, where E is the energy of the heavy quark. Peigne and Peshier (PP) [60] 
have improved this treatment by including the u-channel, which becomes important for 
large energies [see Eq. IA.4l in the Appendix A]. 

For the calculation of radiative energy loss, we consider the treatment of Djordjevic, 
Gyulassy, Levai, and Vitev (DGLV) [611 |57] using opacity expansion, the treatment 
of Armesto, Salgado, and Wiedemann (ASW) [62] using path integral formalism for 
medium-induced gluon radiations off massive quarks, and the treatment of Xiang, Ding, 
Zhou, and Rohrich (XDZR) [63] using light cone path integral approach. Detailed 
expressions for these formalisms are given in the Appendix B. 

5. Results for Energy Loss 

We compare the results for transverse energy loss for a heavy quark using these different 
energy loss treatments for several rapidities. We plot the transverse energy loss of charm 



and bottom quarks, AEt as a function of transverse energy E^ {Jpt + ^^) i^i Figs. [7] 
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Figure 7. CoUisional (dotted lines) and radiative (solid lines) energy loss suffered by 
a charm quark while passing through the QGP 



and [HI at RHIC and LHC energies. 

Several interesting features emerge. We see that the colhsional energy loss for charm 
quarks at RHIC and LHC energies is only marginally dependent on the rapidity and 
the BT formalism gives largest energy loss, as expected. In our treatment, change of 
rapidity implies a change in the temperature of the plasma. Thus these results suggest 
a weaker dependence on the temperature and the average path length for the energy 
loss suffered by charm quarks due to collisions. 

The radiative energy loss, on the other hand, shows a much more complex behaviour 
and is quite different for the different formalisms under consideration. We note that the 
ASW formalism for radiative energy loss gives largest degradation in the energy at all 
rapidities (except for Et < 5 GeV at LHC, where it is comparable to the collisional 
energy loss). We also see that the DGLV and the XDZR formalisms give similar results 
at RHIC energies, at the three rapidities under consideration. On the other hand, at 
LHC energy, the ASW and DGLV formalisms provide nearly identical results for energy 
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Figure 8. Same as Fig. [7] for a bottom quark 



loss for charm quarks at y=0 and y=2, and the corresponding results at y=4 differ by 
about 10%. This, we feel, is due to a more complex dependence on the average path 
length in the ASW formulation. 

The collisional energy loss for bottom quarks using the PP and the Bjorken's 
formulation are seen to be quite similar at RHIC and LHC energies for all rapidities 
under consideration. The BT formulation due to the neglect of the u-channel, gives a 
much smaller energy loss, both for RHIC and LHC energies and at all rapidities. 

We have already mentioned that due to the numerical approximations used, the 
XDZR formulation is not valid for evaluation of the radiative energy loss for bottom 
quarks. The ASW and DGLV radiative energy loss formalisms show a more complex 
dependence on the mass and the average path length. The ASW formulation gives a 
larger energy loss at RHIC energy, though the results are again comparable at LHC 
energy at all rapidities. We note that while the collisional and radiative energy losses 
for bottom quarks at RHIC energy are comparable, the collisional energy loss dominates 
over the radiative energy loss in the Et range under consideration at LHC energy. We 
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have confirmed that at higher Et, the radiative energy loss again starts dominating. 

This rich structure suggests that description of energy loss for one (quark) mass at 
one rapidity, and one energy may not be enough to identify the most reliable treatments, 
for this. 

6. Raa for heavy quarks 

The nuclear modification factor Raa for heavy quarks can be expressed as: 

^^^(^) - TAA{h)da--/d-prdy^ ^'^^ 

where, as mentioned earlier, Taa(&) is the nuclear overlap function for impact parameter 
b, calculated using Glauber model. We get Taa^ 280 fm~^ for Au+Au collisions at 
RHIC and ^ 290 fm'^ for Pb+Pb collisions at LHC, for b = fm. 



u 
P« 



' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' 

Au+Au @ 200A GeV y=() 

y=l 
.-. y=2 




y-" Pb+Pb @ 5.5A TeV 

■■■ y=2 

-■-y=4 




10 15 20 25 30 35 40 
Pt (GeV) 




o 






y=0 Pb+Pb @ 5. 5 A TeV 

y=2 
y=4 




5 10 15 20 25 30 35 40 
Pt (GeV) 



Figure 9. [Left panel] Raa of charm quarks with only nuclear shadowing effect at 
more forward rapidities. [Right panel] Same as left panel for bottom quark. 



As a first step, we give the results for Raa with only nuclear shadowing effect for 
production of charm and bottom quarks at the rapidities considered earlier for RHIC 
and LHC energies (see Fig. [9]). We see that the fairly large masses of the charm and 
bottom quarks, the kinematics, and the rich behaviour of the structure function with x 
and Q^, lends interesting features to Raa- 

We see that Raa for charm and bottom quarks for y = and y = 1 are quite similar 
at RHIC energy. Similarly, the results for y = and y = 2 are only marginally different 
at LHC energy. The results at larger y are more strongly affected due to increased 
variation in the 'x' values (see Eq. |2]) which contribute. In order to do a full justice to 
these interesting results, we now discuss them individually. 

For charm quarks at RHIC energy, we see a suppression at lower p^, an enhancement 
at intermediate pt, and again a suppression at larger pt, for y = and y= 1, while for 
y = 2, Raa starts at about 0.8, goes up to a value slightly more than 1, then drops again 
to about 0.8 at pr~20 GeV, and rise again to beyond 1 at pt~40 GeV. Since the 
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energy loss of the charm quarks always rises with increasing px, this would introduce 
interesting features in Raa after this is accounted for, unless of course the pt spectrum 
for the quarks drops too rapidly. We shall come back to this point again. The increased 
energy at LHC then provides a larger suppression at low pt for all the rapidities. In an 
interesting development, Raa for y = and y = 2 rises beyond 1 at pt~20 GeV, while 
it stays below 1 up to pt ~ 40 GeV, for y = 4. 

The results for bottom quarks are even more interesting. Due to the large mass of 
the bottom quarks, at RHIC energy, the Raa for lower px for y = and y = 1 is already 
starts getting contributions from the region of x where anti-shadowing appears. Thus 
Raa starts at a value which is more than 1 at lower p-r, goes up, up to pr ~ 5 GeV and 
then drops again. For y = 2 on the other hand, it starts at a value close to 1, drops by 
about 10 % at Pt ~ 20 GeV and rises again. At LHC energy, Raa for bottom quarks for 
y = and y = 2 starts at ~30% below 1 and then rises steadily to about 1.1 at px^'^Q 
GeV. The results at y = 4 remain close to 0.8, rising slightly at intermediate px- 
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Figure 10. [Left panel] Raa of charm quarks with the nuclear shadowing effect as 
well as the energy loss at more forward rapidities at RHIC energy. [Right panel] 
Comparison of the relative nuclear suppression using ASW and DGLV formalisms for 
charm quarks at different rapidities at RHIC energy. 



Now let us discuss our results for Raa with the additional inclusion of collisional 
and radiative energy losses. We shall restrict our consideration to inclusion of collisional 
energy loss using the PP formulation and the radiative energy loss using ASW or DGLV 
formulation. 

Fig. [in] gives our findings for charm quarks at RHIC energy. We see that the 
outcome of shadowing and energy loss gives an interesting structure to Raa^ as expected. 
We see that the final Raa starts at about 40 % below 1, goes up to about 0.8 at pr ~ 2 
GeV, and then drops to a value of ~ 0.3 at px ~ 15 GeV. In an interesting development, 
we see that the combination of the shadowing and energy loss gives a marginally larger 
suppression at y = 2 compared to y = 0, even though the fractional energy loss is higher 
at smaller y (see Fig. [7]). We have also given a comparison of Raa by replacing the 
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ASW formulation for the radiative energy loss with the DGLV treatment, and see that 
the former gives a larger suppression at all y (see also Fig. [7]). We recall that the 
single electrons produced from the semi-leptonic decay of charm mesons [26] show Raa 
of about 0.2 ~0.3 for Pt>'^ GeV, in a very encouraging agreement with these results. 
This will be pulbished shortly. 



T 1 1 r- 
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Figure 11. Raa of charm quarks with the nuclear shadowing effect as well as the 
energy loss at more forward rapidities at LHC energy. 



The results for the Raa for charm quarks at LHC energy are shown in Fig. [TTl We 
show the results only for the DGLV formulation for the radiative energy loss, as we have 
seen that it is quite similar to that for the ASW formalism for charm quarks at LHC 
energy. We see roughly similar behaviour, in that the Raa starts from a lower value at 
Pt^Q GeV, rises up to p^ ~ 2 GeV, and then drops to a level of about 0.2 at larger px- 
We also find a marginally larger suppression for y = compared to that for y = 4. The 
results for y = for RHIC energy are also given for a comparison which suggests a much 
larger suppression at LHC, as expected. 

Next we discuss our findings for nuclear suppression for bottom quarks at RHIC 
energy (see Fig. [T2l) . The shadowing and the large mass of the bottom quarks, with 
its consequences, gives an Raa^^-^ at pt~0 GeV, which goes up to about 1.2 at 
Pt^2 GeV, and then drops to about 0.3 ~ 0.4 at larger px- The shadowing results in 
a larger suppression for y = 2 compared to y = (see Fig. [9]), even though the energy 
loss is slightly lower for larger y (see Fig. |8]). Results obtained by replacing the ASW 
formulation with the DGLV treatment show a smaller suppression, as for charm quarks 
(see Fig.[IOD. 

Finally in Fig. [13] we have given our results for Raa for bottom quarks at LHC 
energy for y = 0, 2 and 4 using shadowing, collisional energy loss and radiative energy 
loss using the DGLV treatment. The results using ASW treatment are expected to be 
quite similar as seen from Fig. [HI The results for y = at RHIC energy are also given 
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Figure 12. Same as Fig. [TO]for bottom quarks. 
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Figure 13. Same as Fig. [TT] for bottom quarks. 



for a ready reference. We again see a trend which is common to our results that Raa 
starts at ^0.6 at pr^O GeV, goes up to ~1.1 at pt~2 GeV, and then slowly drops 
to about 0.4 at px ~ 20 GeV. The overall effect of shadowing and energy loss is seen to 
lead to very similar values for Raa from y = to y = 4. Of course the Raa for RHIC 
energy is about twice as large, showing a much reduced suppression. 



7. Summary and Discussion 

We have made a detailed study of charm and bottom production from initial fusion 
of partons in relativistic collision of heavy nuclei. As a first step we have checked the 
usefulness of lowest order pQCD to reproduce the NLO results for the pt distribution of 
charm and bottom quarks in pp collisions. Next we have checked our predictions against 
experimental results of the charm and bottom quarks production from such collisions. 
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We have further obtained the average energy loss suffered by heavy quarks at 
several rapidities for RHIC and LHC energies due to collisions and gluon radiations. 
We have obtained the nuclear suppression factor Raa by additionally incorporating 
nuclear shadowing for the cases under study. A rich picture of dependence of Raa on 
y, pti incident energy, and the mass of the heavy quarks emerges. We have noted that 
our findings would support the suppression of single electrons seen at RHIC. 

Before concluding, we discuss some of the short-comings of the present work. These 
initial calculations can be improved in several ways. We have used (1 + 1) dimensional 
Bjorken hydrodynamics and assumed it to apply at all y. Since the heavy quarks 
will loose most of their energy at very early times, this may not be a serious short- 
coming. Still we are looking at the possibility of using a full fledged (3 + 1) dimensional 
hydrodynamics calculations, also at b 7^ 0. We are incorporating the single electron 
decay of the resulting D and B mesons, along with the results for their back to back 
correlation. We expect this to be rewarding, especially in conjunction with NLO results 
for pp collisions, as it may throw up an interesting detail about differences of NLO 
results and the results with energy loss. These will be published shortly. 

Finally, we conclude that the description for energy loss for one quark mass at 
one rapidity for a particular incident energy may not be sufficient to identify the most 
reliable energy loss treatment for either collisional or radiative energy loss valid for all 
cases. 
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Appendix 

Appendix A. Collisional energy loss 

A.l. Bjorken 

Bjorken argued that the elastic energy loss by partons in the QGP is very similar 
to the energy loss due to ionization due to passage of charged particles in ordinary 
matter [5H]. This treatment was adapted by Braaten and Thoma for heavy quarks |59j . 
The fractional energy loss suffered by the heavy quark due to collisions with the quarks 
and gluons given as: 



dE %TialT'^ ( Nf 



dx 3 V 6 



1 1 - t;^ , 1 + v' 
^-^-log- 

V zv^ 1 — V 



log , (A.l) 

Q'min 



where v is the velocity of the heavy quark. As suggested by Braaten and Thoma we 
use the upper limit of the momentum transfer gmax as yATE and the lower limit of the 
momentum transfer gmin as y/'ivfig. 
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The thermal gluon mass rUg can be expressed as mg=jj,/ \/2 where 



fi=JATT agT'^ (l + -^] is the Debye screening mass. 



A. 2. Braaten and Thoma 

Braaten and Thoma first developed a theoretical formalism to find the collisional energy 
loss of a muon propagating through a plasma of electrons, positrons and photons to 
leading order in QED [M]. This work was further extended by them to calculate the 
collisional energy loss of heavy quarks propagating through QGP |59]. The energy 
loss formulation is given in two energy regimes: The fractional collisional energy loss of 
heavy quarks with energy E « M"^ /T is 
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The fractional collisional energy loss of heavy quarks with energy E » M"^ /T is 
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(A.3) 



where B(f) is a smooth function of velocity having value in the range 0.6-0.7. Braaten 
and Thoma have shown the crossover energy between these energy regimes as -Ecross = 
1.8 X MVT for iV/ = 2. 



A.3. P eigne and Peskier 

In BT formalism, it was assumed that the momentum exchange in the elastic scattering 
process is much less than the energy carried by the heavy quarks. Peigne and Peshier 
pointed out that this assumption is not reliable in the energy regime E » M"^ /T, 
and corrected it in the QED case while calculating the collisional energy loss of a muon 
in QED plasma [65]. This work in QED is then used by them to derive the collisional 
energy loss suffered by heavy quarks while passing through QGP [60] . 

The fractional collisional energy loss suffered by heavy quarks as proposed by Peigne 
and Peshier is 



dE_ 
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andc(A^/) ^ 0.146 iV/ 



0.05. 
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Appendix B. Radiative energy loss 

B.l. Djordjevic, Gyulassy, Levai, and Vitev 

For massless quarks, Gyulassy, Levai and Vitev (GLV) calculated the induced radiation 
to arbitrary order in opacity x" {x = L/X) of the plasma [66] where A is the mean free 
path of the quark. In Djordjevic, Gyulassy, Levai, and Vitev (DGLV) formulation [61], 
the GLV method is generalized to estimate the first order induced radiative energy loss 
including the kinematic effect for heavy quarks. Wicks et al. [57] present a simplified 
form of the DGLV formalism for the average radiative energy loss of heavy quarks: 



where 



1- 



M 



ir A, J I I4E.\ ^ (^2 + ^2)2 



E + p 



(5^ = ml{l - x) + M^x^, 
1 

a: 



Pg^gg + Pg^gg^ 



^g 



'gg 



2/i2 ' 
4 



^qg ~ q ^9gi 

o 1.202 

,, ^o 1.202 



(/32 + K) (/32 Q- +QIQI + Q1 



C = tttW [^'/^' (2^' - ^') + ^' (^' -p')k + Q; (/3^ - 2q'Q; 

+ /.(/3^(-/3^ - 3g^ + /.^) + 2g2g+) + 3/32g2Q,-], 
K = (2px(l - x)f , 

fp = f (/3, q;, q 

U = f (/^, Qt, Q 

and 

f{x, y, z) = i/x4 + 2x2y + z^. (B.2) 
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In Armesto, Salgado, and Wiedemann (ASW), formulation [62] path integral method 
for medium-induced gluon radiation is employed to calculate the radiative energy loss of 
heavy quarks. This formalism provides the analysis of the double differential medium- 
induced gluon distribution by the heavy quarks as a function of transverse momentum. 
The average radiative energy loss is 
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In the above equation the gluon energy, transverse momentum and heavy quark mass 
are expressed as dimensionless parameters. The rescaled dimensionless parameters: 
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We use the parameter no L = 4. 



B.3. Xiang, Ding, Zhou, and Rohrich 

In Xiang, Ding, Zhou, and Rohrich (XDZR) formulation [63], light cone-path integral 
method is used to calculate the gluon radiation from heavy quarks where an analytical 
expression is obtained for the heavy quark radiative energy loss. 
The average radiative energy loss is 
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where a;cr = 2.5 GeV. This analytical expression is derived by using an expansion for 
Bessel function [63j , which is valid only for not too large mass of the quarks. We use 
it only for charm quark. 
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